Solution of a model for the two-channel electronic Mach-Zehnder interferometer 
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We develop the theory of electronic Mach-Zehnder interferometers built from quantum Hall edge states at 
Landau level filling factor v = 2, which have been investigated in a series of recent experiments and theoretical 
studies. We show that a detailed treatment of dephasing and non-equlibrium transport is made possible by 
using bosonization combined with refermionization to study a model in which interactions between electrons 
are short-range. In particular, this approach allows a non-perturbative treatment of electron tunneling at the 
quantum point contacts that act as beam-splitters. We find an exact analytic expression at arbitrary tunneling 
strength for the differential conductance of an interferometer with arms of equal length, and obtain numerically 
exact results for an interferometer with unequal arms. We compare these results with previous perturbative and 
approximate ones, and with observations. 

PACS numbers: 73.23.-b, 73.43.Lp, 03.65.Yz 



I. INTRODUCTION 

Electronic Mach-Zehnder interferometers (MZIs) built 
from quantum Hall edge states have attracted a great deal 
of recent interest, particularly because it has been found that 
they show striking non-equilibrium behaviour>iri2 In these 
mesoscopic devices, quantum Hall (QH) edge states form the 
arms of an interferometer, joined at two points by quantum 
point contacts (QPCs) that serve as beam-splitters. Aharonov- 
Bohm (AB) oscillations in the differential conductance of the 
device are observed 1 as either the applied magnetic field or 
the length of one of the interferometer arms is varied. Ef- 
fects out of equilibrium are probed by studying the visibility 
of AB oscillations as a function of bias voltage: remarkably, 
the visibility does not decrease monotonically with bias, but 
rather shows a sequence of 'lobes' separated by zeros or deep 
minimal 

While theoretical work 14 investigating decoherence in elec- 
tronic MZIs predates these experiments, the problem of un- 
derstanding the origin of the observed lobe pattern has pro- 
vided a fresh focus for such efforts. 15-22 Starting from a de- 
scription of edge states as one-dimensional interacting chiral 
conductors, 23 - 24 two main alternatives have emerged. Onaii is 
that the phenomenon is specific to Landau level filling factor 
v = 2, and arises because in this case there are two types of 
collective modes 24 with different velocities. The otherr^— 
worked out for v = 1, is that the phenomenon is due to mul- 
tiparticle interference effects and requires finite-range, rather 
than contact interactions. 

In this paper we set out an approach that yields exact re- 
sults for a model of an MZI at v = 2 with contact inter- 
actions and arbitrary QPC tunneling probabilities, offering 
a test of previous perturbative 17 and approximate 22 calcula- 
tions. Our approach circumvents a serious technical obstacle 
in the study of QH edge states coupled by QPCs, which is that 
interactions are most easily handled using bosonization,^ 2 ^ 
but this transformation converts the tunneling Hamiltonian 
from a one-body operator in fermionic coordinates to a non- 
linear (cosine) form in bosonic coordinates. Weak tunneling 
(or tunneling probability close to one) can then be treated 



perturbativelyr^i 7 - but additional ideas are necessary in or- 
der to study the situation realised in most experiments, with 
tunneling probabilities close to one half. Two ways around 
this difficulty have been discussed previously. One 21 pro- 
vides numerically exact results, but for simplified models in 
which electrons interact only when they are inside the inter- 
ferometer. The other 22 can be applied to a general model of 
electron interactions, but has so far been implemented within 
an approximation scheme. The complementary route we de- 
scribe here combines bosonization with refermionisation 2 - 7 ^ 
to arrive (for the simplest case of equal-length interferometer 
arms) at a transformed set of fermion creation and annihila- 
tion operators, in terms of which the full MZI Hamiltonian is 
quadratic. Similar methods were used recently by two of us 29 
to treat the related problem 30,31 of equilibration in a QH edge 
state, downstream from a single biased QPC. 

Our main result is that there are important differences be- 
tween the behaviour of a model for an MZI at v = 2 with 
only contact interactions and what has been observed exper- 
imentally. Some of these differences are already apparent in 
behaviour at weak tunneling, computed in Ref. [171 In partic- 
ular for an interferometer with arms of equal length, which is 
the intended situation in most experiments, the visibility of in- 
terference fringes is not suppressed at large bias voltage. One 
might have hoped that the differences would be reduced or re- 
moved outside the weak-tunneling limit, but we show this is 
not the case. We conclude that models with contact interac- 
tions are not sufficient to explain the 'lobe pattern' observed 
in fringe visibility in MZIs. Instead, finite-range interactions 
are required, as discussed in Refs. [l9l - t22l 

The remainder of the paper is organised as follows. In Sec- 
tion [n] we define the model we study for the electronic MZI, 
and in Section [III] we set out the bosonization and refermion- 
isation transformations that we use. We give analytical re- 
sults for interferometers with arms of equal length in Section 
IIVI and describe our approach to interferometers with arms of 
different lengths in SectionjV] presenting numerical results for 
this case in SectionjVI] Finally, in Section lVlIl we discuss our 
results in relation to past theoretical work and experimental 
observations. The numerical procedure used in evaluating the 
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correlators is outlined in Appendix[A] and some comparisons 
with perturbation theory are given in AppendixlBl 

II. THE MODEL 

We consider the model of an electronic Mach-Zehnder in- 
terferometer that is illustrated in Fig.Q] The MZI operates in 
the v — 2 QH regime and is constructed from two quantum 
Hall edges labelled by an index 77 = 1, 2. Each of these edges 
carries two chiral electron channels with spin labels s =\, I. 
Two quantum point contacts a and b, with positions separated 
by distance d\ on the edge 1 and d 2 on edge 2, induce electron 
tunnelling between channels 1J, and 2J,. Two other channels, 
1^ and 2^, are coupled to the rest of the system via contact in- 
teractions. This model was studied using a perturbative treat- 
ment of tunnelling in Ref. and approximately for general 
tunnelling probabilities in Ref. |22[ 

The Hamiltonian H of the model can be written as 

H = H + H tu n, (1) 

where Ho represents the separate edges and H tun describes 
the tunnelling contacts. Defining fermionic fields ^ s {x) that 
create an electron at position x of the channel 77s, with stan- 
dard anticommutation relations 

{^vs( x ):^l's'( x ')} = 5 vv >6 sa >5(x - x'), (2) 
and introducing electron density operators 

p vs (x) = ipl s (x)^ vs (x), (3) 

the edge Hamiltonian (as originally proposed by Wen—) is 

H = -iHvf ^ / tpl s (x)d x ^ vs (x)dx 

T],S J 

n J 

Here, g is the strength of short-range interactions between 
electrons with different spins on the same edge, and short- 
range intra-channel interactions have been absorbed into the 
definition of the Fermi velocity Vf. The tunnelling Hamilto- 
nian H tun = ^tun + ^tun i s characterised by the tunneling 
strengths v a ^ and has the contributions 

i? t Q un = v a e ta ^ u (0)i> 2l (0) + h.c, (5) 
i? t b un = w 6 e^^(d 1 )^4(d a ) + h.c, (6) 

where 4>ab — /3 — a is the Aharonov-Bohm phase due to flux 
enclosed by the MZI. 

A finite bias voltage is modelled by taking sources for dif- 
ferent edge channels to have chemical potentials eV and 
respectively. Two alternative bias schemes have been stud- 
ied experimentally. In the first, which we refer to as single- 
channel bias (SCB), only the channel 14- isjit non-zero chem- 
ical potential, as in experiments of Ref. |2j. In the second, 



which we refer to as two-channel bias (TCB), both channels 
are at the same non-zero chemical potential, as in the 
experiments of Ref. [TH In calculations we establish a steady 
state for the interferometer by using the Hamiltonian H to 
time-evolve from the distant past an initial state \I) that is 
the ground state (or, at finite temperature, thermal equilibrium 
state) of Hq with the specified chemical potentials in each 
channel. 

The edge Hamiltonian Hq has a simple quadratic form after 
bosonization and the difficulties in computing time evolution 
arise from -fftun- For this it is convenient to use the interaction 
representation, taking the tunnelling Hamiltonian H tU n as an 
"interaction" following Ref.[l6|. The time dependence of op- 
erators in this representation, which we distinguish from the 
Schrodinger ones using an explicit time argument, is 

A(t) = e^ ot ie-^ ot . (7) 

Similarly, evolution of the initial state from the distant past to 
time t is induced by the operator 

§(t,-oo) = Texp j-^ H tun (r)dr\ . (8) 

To obtain an expression for the current through the interfer- 
ometer, we recall that the time-dependence of the total density 
on edge 2 is given (for x 7^ 0, e^) by 

'-Y^H^PM] =-(«/ + g)d x £> 2s . (9) 

s s 

Using the continuity equation we can therefore take the cur- 
rent operator to be 

i = e(v f +g)J2p2s(x)\2 (10) 

s 

where x\ < is any point before the MZI, x 2 > d 2 is any 
point after the MZI, and e is the electron charge. This defini- 
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Figure 1 : Schematic view of the model of an electronic MZI treated 
in this paper. The interferometer is constructed from two QH edges, 
1 and 2, each carrying two electron channels, and J,. Horizontal 
arrows indicate the direction of electron motion in the channels. The 
channels 14. and 24. are connected via two QPCs, shown using dashed 
lines. Under single-channel bias, as in the experiments of Ref. |2|, the 
channel 14. is fed from a source at chemical potential /X14, = eV, and 
other channels have sources at p vs = 0. Under two-channel bias, as 
in the experiments of Ref. [TJ sources for the two channels If and 
14- have fj.it — = eV, and the other two channels have sources 
at fi vs — 0. 
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tion is more convenient in our approach that the commonly- 
used alternative in terms of the rate of change of electron num- 
ber on an edge. The steady-state current is then 



I(V) = (I\S\0, -oo)J5(0,-oo)|I) 



(11) 



where the bias voltage V enters through the definition of \I). 
An important characteristic of electronic MZIs which quanti- 
fies the coherence and is measured in experiments, is the vis- 
ibility of Aharonov-Bohm oscillations. It can be expressed in 
terms of the differential conductance 



G(V) = dJ(V)/dV, 



as the ratio 



V(V) 



l (v) + g mia (v)' 



(12) 



(13) 



Here Q max / m i n is the maximum/minimum of the differential 
conductance as a function of AB-phase at a fixed bias voltage. 
Our main aim in this paper is to calculate the visibility at arbi- 
trary bias voltage and QPC tunnelling strengths for the model 
ofEq. ®. 



III. TRANSFORMATIONS 

We use a sequence of transformations to calculate the cur- 
rent in the model for an MZI defined in the previous section. 
Our approach (building on earlier wor k 21,29 by two of us) com- 
bines a treatment of interactions using bosonization with a 
fermionic description of the tunneling at QPCs, and can be 
separated into three steps. 

First, we bring the Hamiltonian Hq into quadratic form 
in a standard way using bosonization. Second, the resulting 
expression is diagonalised by a unitary rotation of bosonic 
fields, where interactions result in appearance of two different 
plasmon velocities. Third, we refermionize the full Hamilto- 
nian H by inroducing new Klein factors. With our choice of 
refermionization transformation, the tunnelling term i? t a un at 
contact a retains its noninteracting form, while H^ un acquires 
an additional phase factor, which can be written in terms of 
electron counting operators. Interestingly, in the case of an 
interferometer with equal arm lengths (d\ = d 2 ) this phase 
factor vanishes, and the interacting problem reduces to a non- 
interacting one, allowing an elementary analytical treatment. 
For an interferometer with unequal arm lengths (di ^ d 2 ) we 
derive an expression for the current in a form suitable for a 
simple and efficient numerical evaluation, and we present the 
details of this approach together with the results. These re- 
sults are for zero temperature, but the expression for current 
through the MZI that we give here are general and can be used 
to study the finite temperature case. 



A. Bosonization 

We start by introducing briefly notation for the bosoniza- 
tion procedure that we use (see Refs. |25[ |26| and HI] for more 



details). We consider initially edges of finite length L, so that 
momentum q is quantized as q = 2nn/L, n E Z, then take L 
to infinity. Fermionic fields if; ns (x) are expressed in terms of 
bosonic fields (j) vs (x) via the bosonization identity 

Vv(x) = (2Tra)- 1 / 2 F vs e l ^^ x e-^ ( - x \ (14) 

where we have introduced Klein factors F vs , particle number 
operators N vs and a short-distance cutoff a. Electron density 
operators are expressed in terms of the bosonic field as 



pr, s (x) = N vs /L- (27r) 1 d x 4>r !S (x). 



(15) 



Bosonic fields have the mode expansion 

ks(x) = -^(27r/ <Z i) 1 / 2 [6 J)s (< Z ) e ^ +h.c.]e-^ 2 , (16) 



9>0 



where bt s (q) is a boson creation operator with momentum q, 
obeying standard commutation relations 

[bvs(q),blf S '(p)} = <W<*««"W (1 7 ) 

From Eqs. (TToT l and (fTTT i the commutation relations for the 
c^rjs (x) operators are 

[$r,s(x),dy4> vs (y)] = -2m8 m >5 ss :8{x- y). (18) 
The mode expansion for the fermonic operator is 

i> vs (x) = ^=J2 d vs(k)e ikx (19) 

and bosonic operators can be expressed in terms of fermions 
via the operator identity 

+00 

% s {q) ^i^jqLfl 2 cl(k + q)c vs (k)- (20) 

k— — 00 

The Klein factors, which add or remove electrons on a given 
channel, have the commutation relations 

{ 77s' f?' 5 '} — 2(5^/ S ss ' , [N^ s , Ffj 1 s 1 ] — ~ ^7777' &ss' -^ris ■ 

(21) 

Using this bosonization prescription, we express the edge- 
state Hamiltonian Hq as 



dx 



2ngh 



27T hv 



E ^Ki + (N ns + 1) . (22) 



'I s 



Similarly we represent the tunnelling Hamiltonian in terms of 
new bosonic fields, obtaining for contact a 



H? un = {2i:ay 1 v a e la Fl i F 2i e l[ ^ 0) -^ m + h.c. (23) 
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and for contact b 

x e *[4>u(di)-4>2i.(d2)} + h c (24) 

Here we have kept the finite contribution to the phase aris- 
ing from the particle number operators and neglected terms 
of higher order in 1/L, which vanish in the thermodynamic 
limit. 



B. Diagonalisation of the Hamiltonian 

The edge-state Hamiltonian Hq given in Eq. d22l i is 
quadratic and can be diagonalized 2 ^ by a unitary rotation 
of the bosonic fields 

(XS+ XA_ XA+ XS-) T = U($if $11 4>2i <p2t) T , (25) 
with the rotation matrix 



1111 

u = -\ 1 - 1 1 - 1 

2 I 1 1-1-1 



(26) 



1-1-1 1 

The transformation preserves bosonic commutation relations 

[Xns( x ), d v x V 's(y)} = -2wiS m ,6 ss rS(x - y). (27) 

Particle number operators N vs transform in the same way as 
bosonic fields, and we use the subscript values rjs = A±, S± 
to distinguish the new operators from the old ones. 
In terms of the new operators, we obtain 



2ir \ hv s » o 2ir . 



where we have introduced the velocities 

v±=v f ± g. 



(28) 



(29) 



The Hamiltonian of Eq. d28t describes four plasmon modes, 
two propagating with velocity u + and two with V- . The final 
term of Eq. d28T i is a constant of motion of H and we omit it 
in the following. 



C. Refermionization 

We refermionize the Hamiltonian H by introducing new 
fermionic fields, expressed in terms of the four new species 
of bosons Xris{x), new Klein factors F vs and particle num- 
ber operators N vs (again employing the subscript values -qs = 



A± , S± to identify the new operators). The bosonization iden- 
tity for the new fields is 

% s (x) = (2na/v s )-^ 2 F ris e i ^^' x e- i ^ x l (30) 
Note that for convenience we normalise to unit current, with 
{4>l s (x),$r,, s ,(x')} =v s 5 m ,5 ss ,8{x-x'). (31) 
The new fermion operators have the mode expansion 

de 



(32) 



where a ns (e) obeys 

{a vs (e),al, s ,(e')}=2Trh5 m ,5 ssl S(e-e'). (33) 
The edge-state Hamiltonian now assumes the simple form 

T]S J T]S J 

(34) 

From Eqns. (Ql and (1341 we obtain an expression for the 
fermionic fields in the interaction representation 

%s(x,t)= f A e -*=(*-x/t..)^( e)< (35) 
J Zirn 

The new Klein factors that enter Eq. ( l30l l can be related to 
the old ones by comparing changes in N ns generated by Klein 
factors in the old and new bases, as described in Refs. |28| and 
l29l An important technical point follows from the fact that a 
unit change in a single N vs in the new basis is equivalent to 
changes of one half in two N vs 's in the old basis; the gluing 
conditions necessary to relate the two Fock spaces have been 
discussed carefully in Ref. |28l A consequence is that only 
bi-linears and not single Klein factors can be transformed be- 
tween bases: the expressions 

HJL=Ht p U> Fs.Fi_=FlF 2t , (36) 
f s- f a+ = F it F n> F s+ F A_ = F it F 2i (37) 

follow directly from the form of the transformation U. We use 
these to transform the combination of Klein factors that enters 
the tunnelling Hamiltonian, finding 



F h F U = FlF lt Fl t F 2l = -F A+ F A _ . (38) 

This gives an expression for the tunnelling Hamiltonian at 
contact a in terms of the new Klein factors 

H? UD = -(2na)~ 1 v a e la Fl + F A _e ll * A + (0) -* A - m + h.c. 

(39) 



7-t 



The crucial property of the transformation of Eq. d261 > is 
that it preserves the form of the tunnelling Hamiltonian at con- 
tact a, in the sense that ±i\A± appears in the exponent with 
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unit coefficient. That makes it possible to use Eq. d30T > to- 
gether with Eq. (T39b to obtain the tunnelling Hamiltonian in 
terms of new fermions, as 



TT<! 

-"tun 



v a e ta ¥ A (0)^ A _(0)+h.c, 



- he^e^^ (di)* A _ (da) + h.c, 



(40) 
(41) 



where we have introduced renormalised tunnelling strengths 

Va,b = -V a ,b/y/V+V-. 

The phase operator Q is given by 



(42) 



Q = j(N2 i d 2 -N n d 1 +N A+ d 1 -N A _d 2 ) 

+ -[4>u{di) - ^2i(d 2 ) - XA + (di) + XA_ (d 2 )] 



Ws+-Afs-]-WA++MA-] 



where 



N vs = / dx p ns (x) 

di 



(43) 



(44) 



Here density operators p V s(x) are defined in terms of the 
bosonic fields Xvs( x ) by analogy with Eq. ( fl"5l ). Eq. d44l 
shows that Af ns counts the number of particles passing the 
QPC a during the time window (—d2/v s , —d\/v a ). 

To find how an applied bias voltage is described in terms 
of the new fermions, note that the chemical potentials in each 
channel enter the grand canonical distribution in the combi- 
nation ^ r)s fi vs N ns . The transformation of N vs therefore im- 
plies that 



(p-s + fj-A-LtA+fj-s-) T = u Vitmmm) 5 



(45) 



which leads to the results shown in Table 1 . It is also interest- 
ing to consider explicitly how an initial state |7) appears in the 
two bases. In the vacuum for the system without tunnelling, 
in which all four channels are at zero chemical potential, all 
single particle eigenstates are filled with electrons from mo- 
mentum — oo up to zero. Klein factors act on this Fermi sea 
by adding or removing particles. Consider as an example the 
case of two-channel bias, for which the initial state has equal 
chemical potentials in both spin channels. To reach this initial 
state from the vacuum we need to add equal numbers of elec- 
trons to each of the channels 1J, and 1^, which is achieved by 
acting repeatedly with the product of Klein factors F^F^. 
Since F^F^ = Fg , this state is also one with equal, 
non-zero particle number in channels A + and S+ and zero 
particle number in the other channels. 

A final step is to re-write the current operator, Eq. ( fTOb . in 
terms of the new operators, as 



I = ev + [p s+ (x) - PA + (x)]\ 



-X2 
-Xl 



(46) 



Equations d34l . ( f40b andl4TTi represent the central result of 
this paper. They give an exact mapping of the initial interact- 
ing problem to the one where the interaction effects have been 



absorbed into the phase shifts for electrons scattering at the 
QPC b. We note that although the channels S± are not cou- 
pled by tunnelling, they generate contributions to the phase 
operators. In a special case of an interferometer with equal 
arm lengths the phase operator Q vanishes and the full Hamil- 
tonian H = Ho + -fftun acquires the form for an MZI without 
interactions but with edges having two different Fermi veloc- 
ities, v+ and v_ . 



IV. INTERFEROMETER WITH EQUAL ARM LENGTHS 

We first discuss the behaviour of the interferometer in the 
special case of equal arm lengths, for which it is possible to 
obtain complete analytical results in a straightforward way. 
With d% = e?2 = d, H tun retains a noninteracting form af- 
ter refermionization. By an elementary calculation, scattering 
between the channels A + and A_ at contact a is described by 
the 5-matrix 



S a 



-ie la t a 

H r 



(47) 



where the transmission and reflection amplitudes are t a = 
sm(v a /h) and r a = cos(v a /ti). The simplification for an 
MZI with equal arm lengths is that scattering at QPC b is rep- 
resented in the same basis by an equivalent matrix Sb- 

To calculate current and hence visibility, consider a fermion 
incident on the interferometer in channel A + and exiting in 
the same channel. Two paths contribute to this process, with 
quantum-mechanical amplitudes A(s\ and Au^ as shown in 
Fig. |2 The amplitudes include distinct phase shifts for each 
path: the phase difference encodes interaction effects and is 
inversely proportional to 

v cS = [i/«_ - 1/u+r 1 . 

We use the form of the current operator given in Eq. d46b . 
At an energy for which particles are incident in the channel 
A + but not in A_, the occupation probability of outgoing 



states in A+ is \A 



- A(jj) | 2 . We hence obtain for the current 



I(V) 



de 
2nh 



{1-\A 



(i) 



A (ii)\ 2 }, 



(48) 



where the integral is taken from the chemical potential in 
channel A_ to that in A + . Defining the tunnelling and reflec- 



tion probabilities T a ^ 
can be written 



ti h and R a , 



the integrand 



1 



R a Tb + T a Rb 

2{T a T b R a R b )? cos(4>ab + ed/hv cS ). 





SCB 


TCB 




+eV/2 


eV 




-eV/2 





fJ-A + 


+eV/2 


eV 




-eV/2 






Table I: Chemical potentials for the new fermions in the two experi- 
mental setups. 
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Aft = r a r b e i£d/Hv - 



A 4 



.4 



A {ll) = -t a t b e^- a h isd f hv - 



Figure 2: Two paths through the interferometer with the correspond- 
ing amplitudes for a particle starting on the channel A+ before the 
QPC a and exiting in the same channel after QPC b. Channels S± 
are not shown. 



Under single-channel bias the energy window for integra- 
tion is from Table I (— eV/2, eV/2). We find for the incoher- 
ent contribution 



e 2 V 

Iq = TrrlRaTb + T a Rb], 

Ann 

and for the contribution sensitive to the phase 4>ab 



(49) 



r (SCB) _ 
1 AB — e 



/eV/2 ^ £ 
—2{T a T b R a R b )i cos(<j> AB +ed/hv cS ). 
-eV/2 



(50) 



From these equations we obtain the conductance 



small tunnelling in Ref. and from an approximate theory in 
Ref. We stress that the result d52b we obtain is in fact an 
exact feature of the model. 

Under two-channel bias the energy window for integration 
is (0, eV) and we obtain for the coherent contribution to the 
current 



t(TCB) 
L AB 



eV 



ds 



and for the conductance 

e 2 



2(T a T b R a R b )i cos((t>AB + ed/hv cS ) 

(54) 



^(TCB) = __ [i?aTb + Tai?h 
Zirn 



2{T a T b R a R b )i cos^ab 



(55) 



In this expression the voltage dependent phase shifts enter the 
oscillating term in a sum with the AB-phase. Because the visi- 
bility is defined as the ratio of the maximum and the minimum 
of conductance at finite bias voltage, the phase <py does not af- 
fect the visibility and the latter is constant and independent of 
voltage 



V( TCB > = V . 



(56) 



To summarize this section, we have obtained exact results 
for the visibility and the phase of Aharonov-Bohm oscilla- 
tions in the conductance of an electronic MZI with equal arm 
lengths at v = 2. The dependence of visibility on bias is 
different according to whether bias voltage is applied to one 
channel or to both channels on the same edge. In the first case 
we obtain oscillations of visibility with a period that is in- 
versely proportional to the arm length; in the second case the 
visibility is bias-independent. In neither case is there a decay 
of visibility with increasing voltage, as found experimentally. 
This demonstrates that a model with only short-range interac- 
tions is insufficient to explain observations. 



g(scB ) = J_ [RaTb + TaRb 

Zirn 

+ 2(T a T b R a R b )^cos(^ v /2)cos<p AB ], (51) 

where <\>y = eVd/hv c s is a voltage dependent phase shift 
which appears as a result of interactions. From Eq. (fSTb we 
find the visibility under single-channel bias 



v (scb) = v |cos(eFd/2fiu e g)|, 



(52) 



where Vo is the visibility of a noninteracting two-paths MZI 
with a single biased channel 



Vo 



2(T a T b R a R b )2 
R a T b + T a R b 



(53) 



We see that, with this bias arrangement, interactions generate 
an oscillatory behavior of the visibility as a function of bias 
voltage with the period inversely proportional to the interfer- 
ometer arm length. This period diverges together with v c g 
in the noninteracting limit. Similar behavior was found for 
an MZI with equal arm lengths using perturbation theory at 



V. MZI WITH UNEQUAL ARM LENGTHS : THEORY 

For an MZI with unequal arm lengths the tunnelling Hamil- 
tonian Eq. ( f40b has a contribution from the phase operator Q, 
and the simple treatment of Sec.[IV]no longer applies. How- 
ever, as we now show, one can still derive an expression for 
the expectation value of the current operator that is amenable 
to a precise numerical evaluation. 

Under the approach outlined in Sec. [II] we require an ex- 
pression for the S-matrix of Eq. ((§). Because 

[K un (n),A 6 un (r 2 )]=0 for n>r 2 

the full S-matrix factorizes into a product of separate S- 
matrices for each contact, with 



where 



S(t, — oo) = S b (t, — oo)S a (t, — oo), (57) 



S a (t, -oo) = exp <{ -- / iJ t a un (r)dr 
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and similarly for §b{t,— oo). Here the usual time-ordering 
can be omitted since H a (r) commutes with itself at differ- 
ent t. This factorization allows us to study the effect of the 
transformations due to QPCs a and b on fermion operators 
separately for each contact. 

Consider first QPC a, described by /J t a un , given after 
refermionization in Eq. (|40T >. We wish to evaluate 

for t = and x > 0. Using the Baker-Campbell-Hausdorff 
formula and writing S a = S a (0, — oo) we obtain 

Sl% A+ (x,0)S a = r a i> A+ (x) - ie ta t a V A _(x) 
Sl$ A _(x,0)S a = r a y A _(x)~ie~ ta t a ii A+ (x).(58) 

Next consider QPC b. To compute the effect of Sb(t, — oo) 
we first use an alternative refermionization scheme, chosen 
so that H b un rather than /J t a un is quadratic. Specifically, we 
introduce a new set of bosons satisfying 

[4>%f(x + di) 4>ii(x + dx) 4>2i{x + d 2 ) (f>2\{x + d 2 )] T 

= U-^+ix) x» A _(x) X b A+ (x) X b s .(x)} T . (59) 

We then define new fermions related to x^ s bosons in the 

same way as in Eq. ( f30b . Now H b un expressed in terms of b 
fermions is 

H b UD = v b e if) i> A \ (0)$ b A _ (0) + h.c, (60) 

and so the transformations for the $ A± operators due to con- 
tact b have the same form as Eq. d58l ). The two sets of bosonic 
fields Xns (x) and X b ( x ) 316 related by 

X h n ±{x) = )^[Xr,±{x + di) + Xr-i±(x + d 2 )} 

+ l;[Xv±(x + d^ - Xv±(x + d 2 )}, (61) 

where ?; = A and rj = S or vice-versa. 

The current operator in terms of the 6-densities is simply 

i = ev + \p b s+ {x)-p\ + {x)]\ll^l (62) 

and the next step is to find the transformation induced by 
§b(0, — oo) = Sb on the quantities that enter the current oper- 
ator. To do this, we express the densities in terms of normal- 
ordered fermion fields in the standard fashion, as 

v s p TIS {x) = :i>l s {x)* r , s {x)t, (63) 

where normal ordering * ... * is defined with respect to the vac- 
uum state |0), which obeys 

a vs (u)\Q) = 0, lu > 0, 
o+ (w)|0) = 0, w<0. 



The transformation can then be derived using the equivalent 
of Eq. 08] ) for b fermions, and we obtain for x > 

Slv + [p h s+ (x) - p b A+ (x)]S b 

= v+p b s+ (x) - [r 2 b v+p b A+ (x) + t\v-p\_ (x)\ 

+ \ie^rbib^\{x)¥> A _{x) +h.c], (64) 

while for x < 

Stlp b s+ (x) - p b A+ (x)]S b = {p b s+ (x) - p b A+ (x)}. (65) 

The terms on the middle line of Eq. d64l i give the incoherent 
(or AB phase-independent) contribution Iq to the steady-state 
current I(V) while those in the final line make the coherent 
contribution I A b- 

To complete the evaluation of the steady-state current, we 
return to the operators i& vs (x) using the transformation of 
Eq. doTt . In this way we find 

I = ev + (I\p A+ (xi)\I) 

- e(I\St[rlv+p A+ (x2)+tlv-p A _(x 2 )]S a \I}. (66) 

Similarly we obtain for the coherent contribution 

I AB =e[ie^r b tb{I\Sl^^{d 1 )^ A _{d 2 )^I)+h.c.}. 

(67) 

These constitute the required expressions for the current 
through the MZI. As we show in the following sections, they 
are suitable for numerical evaluation. They would also pro- 
vide the starting point for an approximate analytical treatment, 
although we do not explore that direction here. 

VI. BEHAVIOUR WITH UNEQUAL ARM LENGTHS 

There are two aspects to the behaviour of an MZI with con- 
tact interactions and equal arm lengths that are strikingly dif- 
ferent from what one might expect in more general models: 
first, there is no suppression of the visibility of interference 
fringes in the differential conductance at high bias voltage; 
and second, with two-channel bias, visibility is completely 
independent of bias. Both these aspects change when arm 
lengths are unequal. In this section we discuss the physical 
reasons for these changes and present detailed numerical re- 
sults. 



A. Qualitative discussion 

Suppression of interference fringe visibility at high bias 
for an MZI with d\ ^ d 2 is due in our treatment to fluc- 
tuations in the phase Q, appearing in Eq. d67b . The bias 
dependence of these fluctuations arises via the contributions 
from J\f A+ and Af A _ to Q [see Eq. d43ll. In detail (taking 

for definiteness d 2 > d\ and p A+ > p A _), Af A + counts 
particles in the channel A + at time zero that passed QPC a 
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in the time interval (— d,2/v+, — d\/v+) without tunnelling, 
while Ma- counts those that did tunnel in the separate in- 
terval {—d%/v-, —dx/v-). Bias dependent fluctuations in Q 
come from particles that pass QPC a during the portions of 
these time intervals that do not overlap, since the contribution 
of such particles to Q depends on whether they tunnel. De- 
veloping this picture, one can identify two separate regimes, 
according to the value of the ratio 7 = d\v+ jd^v-. For 7 < 1 
it is the intervals (—d\/v-,—d\/v+) and {—d<2,/v-,—d<2,/vj r ) 
that contribute to fluctuations of Q, while in the oppo- 
site case 7 > 1 it is the intervals (— d\/v+, — ^2/^+) and 
(—d\/t>-, —dz/v-). The scale for the bias voltage at which 
interference is suppressed is the one at which the fluctuations 
in particle number within an energy window eV and on the 
given spatial intervals are 0(1). If T a ~ 1/2, this scale is 
V ~ e±/e, where for 7 > 1 the relevant energy is 



e+ = H{(di + da) [l/«_ - !/«+]}" 



and for 7 < 1 it is 

e_ = h{(d 2 -di)[l/«- 



!/«+]}■ 



(68) 



(69) 



As expected, this voltage scale diverges both in the non- 
interacting limit and for equal arm lengths. 

A related argument can be used to understand why, for two- 
channel bias, oscillations in visibility occur only with arms of 
unequal length. In this case the relevant feature is the bias- 
dependence of the average value of Q, rather than its fluctua- 
tions. For two-channel bias we find 

pT V 

n(I\StQS a \I) = (d 1 -d 2 )—^V = -. (70) 

The effect of this bias-dependent phase can be modelled by 
including it in the integrand of Eq. i54i . making the replace- 
ment <))ab — > 4>ab + V/Vq. Then for eVo ^> hv e s/d the 
visibility is 




2 3 

eVdi/2jrSv e H 

Figure 3: Single-channel bias: dependence on arm length difference 
and interaction strength. Top: Normalized visibility as a function of 
bias voltage for T a = Tt = 1/2. Black dashed line (fe/di = 1, 
g/vf = 0.25; red solid line dijd\ = 1.1, g/vf = 0.25; blue dashed 
lined2/di = 1.2, g/vf = 0.25, green dot-dashed line tfe/di = 1.2, 
g/vf = 0.125. Bottom: Phase of AB-fringes as a function of bias 
voltage for the same parameters. 



V~ V 



2hv cS 

1 — — cos 4> v 

eVod 



(71) 



It has oscillations as a function of bias, with an amplitude that 
vanishes for equal arm lengths. 



B. Numerical results 

In this section we present numerically exact results for the 
visibility and phase of interference fringes in the differential 
conductance of a MZI, obtained by evaluating Eq. d67l i using 
the methods outlined in AppendixlAl 



c?2 = d\ the visibility oscillates with constant amplitude as 
a function of bias voltage, following Eq. d52l . and the phase 
of AB-fringes shows exact steps of height tt at the zeros of 
visibility. The energy scale of oscillations is Ed = Hv e g jd and 
diverges in the non-interacting limit. For d\ 7^ d2 the visi- 
bility develops a decaying envelope on the energy scale e + or 
e_ (depending on the value of 7), minima in visibility are no 
longer exact zeros, and the steps in the phase as a function of 
bias are smooth. 



1. Interferometer with single channel bias 

The dependence of interference on difference in arm 
lengths and on interaction strength is shown for an MZI with 
single-channel bias and T a = T\, = 1/2 in Fig. [3] For 



Variations with the tunnelling probability T a are shown 
in Fig. |U When the QPC a is tuned away from the half- 
transparency, the visibility and the phase of AB-fringes be- 
come irregular functions of bias voltage. Changes in X& alter 
only the overall scale for the visibility, and not the form of its 
dependence on bias. 
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eVd 1 /27rAv e fj 



Figure 4: Single-channel bias: dependence on T a . Top: Normalized 
visibility as a function of bias voltage for d%/d\ = 1.2 and g/vj — 
0.25, T b = 0.5 and different T a . Red solid line T a = sin 2 tt/8; blue 
dashed line T a = sin 2 37r/8; green dot-dashed line T a — sin 2 ty/6. 
Bottom: Phase of AB-fringes as a function of bias voltage for the 
same parameters. 



2. Inteferometer with two-channel bias 

The dependence of interference on difference in arm 
lengths and on interaction strength is shown for an MZI with 
two-channel bias and T a = T b = 1/2 in Fig. [5] For equal arm 
lengths the visibility is independent of bias and the phase of 
AB-fringes varies linearly with voltage, following Eq. d55l ). 
For di ^ d,2 the visibility develops small amplitude oscilla- 
tions and is suppressed at large bias: the voltage period of 
oscillations is consistent with the value 2irhv c s/ed expected 
from Eq. ( |7TI ). The variation of phase with voltage is no longer 
exactly linear but shows no well-defined steps. 

Results for several values of T a are displayed in Fig. [6] os- 
cillations of visibility with bias are irregular but some well- 
defined minima develop for values far from T a = 1/2. 

3. Suppression of visibility at high bias 

From the discussion in Section lVI Al we expect the voltage 
scale for the suppression of visibility at high bias to be set by 



1.4 - 



1.2 - ,'\ 




eVdi/27rftv e ff 



7 - 



6 - 




eVdj/^TT/iVrff 



Figure 5: Two-channel bias: dependence on arm length difference 
and interaction strength. Top: Normalized visibility as a function of 
bias voltage for T a — Tb = 1/2. Black dashed line d^jd\ = 1, 
g/vf — 0.25; dot-dashed violet line d-2,jd\ = 1.05, g/vf = 0.25; 
red solid line dz/d\ = 1.1, g/vf = 0.25; blue dashed line d2/di = 
1.2, g/vf = 0.125. Bottom: Phase of AB-fringes as a function of 
bias voltage for the same parameters. 

e±. Behaviour consistent with this is shown in Fig. [7] here 
7 > 1 for all parameter sets, and £+/e sets a common voltage 
scale to the envelope for visibility oscillations. We note that 
for T a — 1/2 this envelope is approximately Gaussian. 



VII. DISCUSSION 

In this section we discuss our results in relation to experi- 
mental observations and previous theoretical work. To sum- 
marise briefly: building on techniques developed in earlier 
studies 21 '^ 9 we have presented exact analytical and numeri- 
cal results for the non-equilibrium behaviour of an electronic 
Mach-Zehender interferometer built from quantum Hall edge 
states at filling factor v — 2, using a model with only contact 
interactions. A key feature of the results is that, for an interfer- 
ometer with nearly equal arm lengths, two scales are present 
in the dependence of fringe visibility on bias voltage: oscil- 
lations in visibility have the voltage period 2irv e g/ed, with 
d « di w g?2, while their envelope falls off on the larger scale 
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1.4 - 




eVdi/27rAv e ff 



Figure 6: Two-channel bias: dependence on T a . Top: Normalized 
visibility as a function of bias voltage for (fe/di = 1.2, g/vf = 
0.25, Tb = 1/2 and different T a . Red solid line T a = sin 2 tt/8; blue 
dashed line T a = sin 2 37r/8; green dot-dashed line T a = sin 2 7r/6. 
Bottom: Phase of AB-fringes as a function of bias voltage for the 
same parameters. 



2wv e s /e\d\ — d% \ . A further feature is that oscillations in vis- 
ibility are much clearer in one bias scheme (single-channel 
bias) than in an alternative scheme (two-channel bias). Fi- 
nally, AB fringes show a varying degree of phase rigidity un- 
der small changes in bias: for single-channel bias it is per- 
fect when arm lengths are equal, but decreases with increasing 
\d\ — d%\\ for two-channel bias it is absent if arm lengths are 
equal, and otherwise at most limited. 

Some but not all aspects of these results match observa- 
tions. Most importantly, the originally reported 2 iobe pattern' 
is reproduced by calculations for the experimentally appropri- 
ate single-channel bias scheme with d\ ^ d2 (see Fig. [3]). In 
addition, as previously discussed in the perturbative context, 17 
the differences in behaviour between the two bias schemes 
(compare Ref. Q with Ref. [l2h are matched by differences in 
calculated behaviour (compare Fig. [3] with Fig. |5j. On the 
other hand, as at most a few oscillations of visibility are ob- 
served in the lobe pattern^^ the voltage scales for oscillations 
and for their envelope are not well-separated: since the sam- 
ples concerned are intended to have almost equal arm lengths, 
this is at variance with the behaviour of the model we study. 



Moreover, intentional changes in the length of an arm appear 
to have a much smaller effect experimentally 2 than in our cal- 
culations. 

It is interesting to go beyond these qualitative comparisons 
and attempt an estimate of the key theoretical parameter, v e s, 
which characterises interaction strength in our model. This is 
possible using measurements by Roulleau and collaborators, 
described in Ref. 9. The relevant experiment, in the notation 
of our Fig.[T] involves applying separate biases to the channels 
1 f and 1 4_, with the other two channels grounded. Specifi- 
cally, applying a voltage to 1 f, this channel acts as a modu- 
lation gate, changing the phase of AB oscillations in conduc- 
tance. From Fig. 5. 17a of Ref.lKJ a bias of 49/iV generates a 
phase shift of 2tt in a sample with d% = d-2 = 11.3/im. Us- 
ing the results described in Section [TV] we obtain from these 
data the value v e g = 6.7 x 10 4 ms _1 . Remarkably, this is 
very close to the estimate (v e g = 6.5 x 10 4 ms -1 ) obtained 
in Ref. 2^ from a theoretical fit of an experiment on equilibra- 
tion of QHE edge states, although in general u e g is expected 
to vary with sample design and magnetic field strength. 

Our results should also be compared with a body of earlier 
theoretical work, which (apart from the early study of Ref.fl5h 
can be separated into investigations of the effects of long range 
interactionsr^r— and calculations for the mode l 17 ' 22 with 
contact interactions that we have studied here. The initial 
treatment of this modeU^ was perturbative in the tunnelling 
amplitude at the two QPCs, and so appropriate for T a j <C 1 
or 1 — T a {, <C 1, while an approximate non-perturbative ap- 
proach has been described in Ref. |22|. The advance we have 
presented here is to handle tunnelling exactly. The successes 
and weaknesses of the perturbative calculations 17 in account- 
ing for observations are similar to the ones we have described 
above. In particular, the fact that the voltage scale for suppres- 
sion of AB oscillations is set by the difference in arm lengths, 
and diverges for d-y — c?2 appears at the perturbative level. 
Knowing only the original, leading order result, 17 one might 
have hoped that contributions higher order in T a would elim- 



o. 




0.05 0.1 0.15 0.2 0.25 



(cV(d 2 - d 1 )/2^fiv dr ) 2 

Figure 7: Suppression of coherence at high bias: visibility on a log- 
arithmic scale as a function of bias voltage, two-channel bias: red 
solid line dz/di = 1.1, blue dashed line dz/di = 1.2. Single- 
channel bias: green dotted line dz/di = 1.1, violet dot-dashed line 
d 2 /di = 1.2. In all cases g/vf = 0.25 and T a = T b = 1/2. 
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inate AB oscillations at large bias for all d\ — d% (as happens 
when shot noise is introduced using a separate QPC 18 ). The 
calculations we have described show (in agreement with the 
approximation of Ref . 1221) that this is not the case. 

An implication of the work we have presented is that one 
must take account of long range interactions in order to under- 
stand in full experiments on non-equilibrium effects in MZIs. 
Calculations that do include long range interactions^— suc- 
cessfully reproduce many aspects of the observations, but 
have been done for v = 1 while much of the published 
data are for v = 2. Approaches such as the one of Ref. |22j 
that include both long range interactions and the two chan- 
nels present at v — 2 are therefore desirable. In this context, 
the results we have presented provide a testing ground for ap- 
proximation schemes. Interestingly, the same model for QH 
edge states that we have studied here, of contact interactions at 
v = 2, seems to account more successfully for experiments on 
relaxation of a non-equilibrium electron distribution 29 ' 30 than 
for the behaviour of an MZI with equal length arms. While we 
have no detailed understanding of this difference, we point out 
that AB interference is a more sensitive probe of the QH edge 
than a measurement of the electron distribution. In particular, 
if as is likely, weak long range interactions are present in addi- 
tion to contact interactions, they can be expected to suppress 
interference in a MZI at high bias 16 without much changing 31 
the observed relaxation process. In outlook, we hope that it 
will be possible to extend the techniques we have developed 
here to treat other phenomena in edge states far from equilib- 
rium. 
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Appendix A: Numerical evaluation of the tunnelling 
conductance 

Here we outline the numerical procedure we use to obtain 
results for an interferometer with unequal arm lengths. The 
task is to evaluate Eq. d67l i. in which the average is taken with 
respect to the scattering state produced by QPC a. This scat- 
tering state is generated by the action of the operator S a on a 
Slater determinant describing filled Fermi seas that are defined 
by chemical potentials A±. The scattering operators &l . . . S a 
appearing in Eq. ( |67| | can more alternatively be taken to trans- 
form the operators they enclose rather than the state vectors, 
and the averages we require then have the general form 

Aij = (ate^w^'^'aj), (Al) 

where a\ is an electron creation operator and the index i labels 

both channel and energy eigenstate. Since Y^kl ^kloX&l is 
quadratic in fermion operators, we can use Wick's theorem to 



obtain 

= detM x[hx M'^ji, (A2) 

with 

M = I + (e l " - I)n . (A3) 

Here I is an identity matrix and h is the density operator, 
which can be written in the basis of scattering states 

Uij = (i\n\j) = (otoi). (A4) 

The matrix M has a similar form to those appearing in the 
problem of full counting statistics 33 and in nonequilibrium 
bosonization. 34 We evaluate Eq. (IA21 > numerically using an 
energy eigenstate basis for a system of finite length L with 
the eigenstates on a given channel occupied from the bottom 
of an energy window up to the corresponding chemical poten- 
tial. We use at most about 2000 states in each channel. 

As a check we present in Fig. [8] a comparison of results 
from these numerical calculations at weak tunnelling with 
those from perturbation theory in tunnelling strength. The 
agreement is essentially perfect on the scale visible in our fig- 
ures. 



Appendix B: Perturbation theory in the small tunnelling limit 

In this appendix we recall the results of perturbation theory 
in tunnelling strength and compare behaviour for T a ^ = 1/2 
with that at weak tunnelling. 

Perturbation theory in tunnelling strength was applied to 
MZIs at v = 1 in Refs.QJ and^and at v = 2 in Ref.[H 
For the model we study, the current through the interferometer 
at small v a , Vb is given by 

1{ y) = J 1 [{v « +v 2 b )e- ieVt/H iM9(0,t) 2 } 

+ {v a v b e^ AB e~ ieV ^-^/ h ilm[g(d l7 t)g(d 2 ,t)} 

+ c.c.}], (Bl) 




Figure 8: Comparison of the exact result (circles) with perturbation 
theory (solid line) in the small tunnelling limit for g/vf = 0.75 and 
d 2 /d! = 1.2. 
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where g(x, t) is 



eVdi /2irhv e $ 

Figure 9: Comparison of the exact results at T a = Tt = 1/2 (sym- 
bols) with those at weak tunnelling (lines). Top: single-channel 
bias; bottom: two-channel bias. Solid black line and red circles: 
d-zjd\ = 1.1; dashed black line and blue triangles: (fe/di = 1.2. 
All data for g/v f = 0.25. 



g(x,t) 



2ir (x - v + t + ia) 1 / 2 (x - v_t + ia) 1 / 2 



Hereto = div/v+v~ for single-channel bias, and to = d±/v+ 
for two-channel bias. The integral in the Eq. (IB lb can be 
evaluated 17 analytically in the strong coupling limit v+ /v- ^> 
1, giving (for example) with two-channel bias the visibility 



V, 



(TCB) 
port 



(V) =V \Jo(eV[d 2 - d!]/2hv-)\, 



(B2) 



where Jq(x) is a Bessel function. For general interaction 
strength and arm lengths one can evaluate Eq. (IB 1I > numeri- 
cally. We compare the results of these calculations with our 
results at T a = T& = 1/2 in Fig. |9l using the same param- 
eters as in Figs. [3] and [5] We note that while the qualitative 
behaviour is similar for both strong and weak tunnelling, spe- 
cific features differ greatly, especially at large bias. These dif- 
ferences appear to be much greater than those found in the ap- 
proximate treatment of finite tunnelling strength 22 (compare 
ourFig.[3]with Fig. 5 of Ref.[H 
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